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ABSTRACT:  Modeling  of  a  flexible  robot  arm  mounted  at  a  prismatic  joint  leads 

to  a  moving  boundary  problem.  In  a  control-theoretic  context  with  joint  motion  as 
control,  these  become  free  boundary  problems.  Some  problems  are  discussed  for  trans¬ 
verse  vibration  (the  beam  equation)  but  the  principal  result  is  an  exact  controllability 
theorem  for  longitudinal  vibration  (the  wave  equation). 

Key  words:  free  boundary,  robot  arm,  prismatic  joint,  beam  equation,  feedback,  wave 
equation,  exact  controllability 


1.  INTRODUCTION 


We  consider  the  modeling  of  a  robot  arm  consisting  of  a  flexible  beam  with  a  prismatic 
joint  —  that  is,  a  mounting  of  the  arm  in  which  one  end  is  clamped  but  permits 
controlled  in-out  motion  so  the  extension  of  the  arm  (i.e.,  the  relevant  portion  of 
the  beam)  has  variable  length  f  —  f.(t).  In  this  case  ‘flexible’  means  only  that  the 
dynamics  of  the  beam  must  be  modeled  by  a  partial  differential  equation:  treatment 
as  a  rigid  body  is  an  inadequate  approximation  in  that  vibration  (i.e.,  time-dependent 
deviation  from  the  nominal  rigid-body  position)  is  a  significant  aspect  of  the  problem. 
The  space-time  domain  on  which  we  consider  the  dynamics  is  thus 


-  ! It  : 


{(/,.s)  :  0  <  t  <  7’;  0  <  ,s  <  C(t)} 


'This  research  has  been  paiti-illv  c«tppoi(ed  by  the  H.S.  Air  Force  Office  of  Scientific  Research 
under  grant  #AFOSR-87  moo 

2This  material  was  presented  at  the  international  colloquium:  Free  Boundary  Problems:  Theory 
and  Applications  held  at  Irsee  (Havana,  Germany)  June  11-20,  1087.  The  present  text  will  appear  in 
the  proceedings  of  that  conference,  edited  by  K.-ll.  Hoffmann  and  J.  Sprekels. 


where  s  denotes  length  (in  material  coordinates,  measured  from  the  tip  of  the  arm). 
This  is  certainly  a  moving  boundary  problem.  From  a  control-theoretic  point  of  view, 
we  note  that  the  significant  question  is  precisely  the  determination  of  the  control  £(•) 
so  we  will  have  a  free  boundary  problem. 

In  the  next  section  we  analyze  purely  longitudinal  vibration  of  such  a  beam.  This 
is  not,  of  course,  the  problem  of  greatest  practical  concern  but  we  will  be  able  to 
treat  it  with  some  degree  of  completeness  —  showing  exact  controllability,  subject 
only  to  some  rather  mild  conditions.  The  more  significant  problem  of  treating  trans¬ 
verse  vibration  in  this  context  requires  more  extensive  analysis  than  can  be  presented 
here.  A  brief  discussion  of  such  problems3  is  presented  in  Section  3.  For  our  present 
purposes  we  take  the  beam  to  be  (nominally)  horizontal  along  the  x-axis  of  the  ‘lab¬ 
oratory  coordinate  system1  with  the  joint  rigidly  fixed4  at  the  origin.  Included  in  our 
model  is  a  load  of  mass  m  supported  at  the  free  end  of  the  arm.  We  assume  that  the 
system  is  (except  for  the  boundary  condition  at  the  joint)  conservative  —  dissipative 
mechanisms  such  as  friction  and,  e.g.,  air  resistance  can  be  neglected. 


2.  EXACT  CONTROLLABILITY 

Here  we  restrict  our  consideration  to  the  purely  longitudinal  vibrations  of  the  arm. 
Suppose  we  let  x(t,s)  be  the  actual  position  (in  the  laboratory  coordinate  .system) 
of  the  material  point  s  at  time  t.  The  formulation  of  the  dynamics  takes  a  much 
simpler  form  if  we  take  the  unknown  to  be  v  (s  I  x]  so,  e.g.,  v  —  t  at  the  joint.  For 
a  hyperelastic  stress-strain  law6  one  can  easily  show,  using  Hamilton's  Principle  of 
St  ationary  Action  and  a  bit  of  manipulation,  that  v  must  satisfy  the  one-dimensional 
wave  equation: 

(2.1)  pvlt  -  !/»(•,«,),  on  Qt 
with  the  boundary  conditions: 

(2.2)  v  -  f  at  the  joint:  -s  =  i  =  £((.), 

(2.3)  rn  v,,  ~  at  the  tip:  s  —  0. 

■’Work  on  such  problems  has  Itch  initiated  l>y  P.K.O.  W.iiir  and,  indeed,  |2|  should  be  acknowledged 
as  the  principal  stimulus  to  the  present  work.  Many  interesting  and  important  problems  remain  open. 

,,ln  practice  one  would  expect  that  the  joint  may,  itself,  be  mounted  on  a  moving  support  (e.g., 
on  another  arm)  and/or  that  theie  i«  coupling  with  a  controlled  rotation  around  a  vertical  axis.  One 
might  also  wish  to  include  the  po««il>ililv  that  the  mounting  is  nominally  rigid  yet.  subject  to  vibration 
giving  'noise’  in  the  boundary  condition  theie;  compare  |l|. 

r’ We  make  the  usual  assumption  that,  the  elastic  potential  ♦  —  ♦(.*,  r)  is  convex  in  r  with  minimum 
at  ^(.*,0)  =0,  We  set  ^(s,  r)  -  't'.(vr)  and  let  p(.«)  be  the  mass  density. 


To  this  we  adjoin  specification  of  initial  conditions: 


(2.4)  £(0)  =  £°  and  v  =  v°,  vt  =  w°  on  [0,£°]  at  t  =  0, 

subject  to  (2.2):  v°(£°)  =  £°  as  a  compatibility  condition. 

Suppose,  in  addition  to  (2.4),  one  were  also  to  specify  a  target  state  for  the 
terminal  time  t  —  T: 

(2.5)  C(T)  =  CT  and  v  =  vT ,  vt  =  wT  on  [0,  tT\  at  t  =  T, 

again  subject  to  (2.2).  Our  object  is  to  prove  exact  controllability: 

There  should  be  a  time  T  >  0  and  an  admissible  control  function  £(•)  on 
[0,r|  such  that  (2.1),  (2.4),  (2.2),  (2.3)  give  (2.5) 

under  suitable  (mild)  hypotheses.  By  “admissible1’  we  mean  that  £(0)  =  £°,  t[T)  =  £r, 
and0 

(2.6)  (€(£,)  -  £(/2)|  <  k|<,  -  l2\  with  k  <  I 

which  will  ensure  that  the  free  boundary  curve:  F  :=  {s  =  £(<)}  is  ‘time-like’  so  a 
characteristic  |.s  ±  t  =  const.]  can  intersect  F  only  once  (transversely). 

For  simplicity  of  exposition  we  present  in  detail  only  the  linear  case7  for  which 
(2.1)  reduces  to 

(2.7)  vtt  =  t;„  on  QT- 

In  this  case  the  method  of  characteristics  gives  the  simple  decomposition 

(2.8)  v[t,s)  =  f{s  +  t)  +  g(s  -  l) 

into  left-  and  right-moving  waves.  Note  that  differentiating  (2.8)  gives  the  identities: 

(2.9)  f'(s  +  t )  =  (t>,(M)  +  u((«,.s)|/2, 

g'(.s  -t)  =  (v,(/., .s)  -  v,(t,s)]/2. 

The  representation  (2.8)  gives  the  general  (weak)  solution  of  (2.7)  and,  in  terms 
<>r  /,  <7,  the  boundary  conditions  (2.2),  (2.3)  now  become: 

(2.10)  T<)  +g{l{l)  -t)  =  £(/.), 

(2.1!)  m\f"{s)  t-  g"{  -s)]  =  f'{s)  4  g'{- s). 

0  Physically,  ill ie  corresponds  to  a  requirement  that  the  controlled  motion  at  the  joint  must  always 
be  slower  than  the  wave  propagation  speed  along  the  beam.  For  smooth  t(  )  it  is  equivalent  to 

!*(«)!<*<  i. 

7l.e.,  'J»(s,r)  =  ‘/9(s)r2  so  ip  =  dr.  We  further  simplify  by  assuming  the  beam  is  homogeneous  so, 
in  suitable  units,  p  =  0  =  1. 


Note  that  we  can  differentiate  (2.10)  with  respect  to  t  to  obtain  an  ODE  for  £: 


(2.12) 


;  /'(<  +  t)  -  „'(<  -  I) 

"  i  -i/’(< +<)-?’(<  -or 


From  the  inequality.  |u  -  v|  +  |u  4-  t;j  <  2  max{|u|,  |t>|}  it  would  follow  that  (2.12)  is 
well-behaved  —  indeed,  giving  (2.6)  —  if  we  were  always  to  have 


(2.1.1) 


where  defined. 

The  initial  data  v°,w°  in  (2.4)  determine  /, g  on  (0,  £°]  using  (2.9)  with  <  =  0.  In 
particular,  we  have  g'  determined  on  Dq  :=  {(t,s)  :  0  <  s  -  t  <  £0}.  We  then  can  use 
(2.11)  as  an  ODE  for  g'  to  obtain: 

(2-14)  g'(s)  =  +  w°|(-s)  +  e^mw°{0) 

+  ^f  *  e(,+r>/m(v?  +  w°\{r)dr 

for  -£°  <  s  <  0,  which  defines  g'  from  (2.4)  on  D\  :=  ((<,a)  :  —  £°  <  s  -  t  <  0}. 
Finally,  we  note  that  w(£°,0)  is  determined  by  the  data  of  (2.4). 

Similarly,  we  note  that  the  target  data  in  (2.5)  directly  determine  /'  on 
{(<,«)  :  T  <  s  +  t  <  tT  T  T},  using  (2.9)  with  t  =  T,  and  give,8  now  using  (2.11)  as 
an  ODE  for 


(2.15)  f{T  -  r)  =  l|„r  +  »r)(r)  + 


for  0  <  r  <  £T,  which  defines  /'  from  (2.5)  on  V{  : =  {(<,$)  :  T  -  lT  <  s  +  t  <  T }. 
From  (2.5)  we  also,  then,  determine  v(T  -  £T,0). 

We  now  choose  Cauchy  data  on  the  segment:  tj  :  —  {(CO)  :  £°  <  t  <  T  —  £T}, 
taking  v(<,0)  linear  in  t ,  making  v,  constant: 


(2.16) 


v,(C0)  =  c  := 


v(r  -  £T,0)  -  r(£n,0) 


T  -  |£°  +  £r] 

on  §.  This  gives  v,(t,0)  —  tnvtt  =  0  by  (2.3)  with  t/>(r)  =  r.  Using  (2.9),  we  now 
have 

(2.17)  /'  =  c/2  on  D{  := {(t,s)  :  £°  <  .?  +  /  <  T  —  £r}, 

g'  =  -c/2  or,  V{  :=  {(<,.?)  :  lT  -  T  <  9  -  t  <  -£"}. 

"To  avoi<l  any  possibility  of  conflict  we  assume  t  hat  T  is  chosen  so  T  >  ±  tT ,  giving  <  (T  -  T 

and  ('  <  T  +  tT ;  later  we  will  strengihen  this  requirement  to  (2.18).  Note  that  the  computation,  here, 
from  (2.5)  is  independent  (to  within  translation  in  t)  of  the  particular  choice  of  T. 


Note  that  if  we  require  that 


(2.18)  T  >  (£°  +  lT\  4-  |t/(T  -  t\0)  -  i/(£°,0)|/k, 

then  (2.16),  (2.17)  give  (2.1.1)  on  P{,P{.  If  we  now  impose  the  requirement  that 

(2.19)  k?|>°|  <  k/4  on  (0,£°), 

kTi>ri  < °n  mti, 

then  (2.13)  holds  on  ,  Pi  directly  from  (2.9)  and  it  easily  follows  from  (2.14),  (2.15), 
that  we  will  also  have  (2.13)  on  P{,Pf. 

Theorem  1:  Suppose  one  is  given  the  data  of  (2.4),  (2.5)  subject  to  compatibility 
with,  (2.2)  and  the  condition  (2.19)  with  k  <  1.  (Assume  w°,  ivT  are  continuous 
and  v'\  vT  G  C1.)  Then  there  exists  T  >  0  and  an  admissible  control  function 
£(•)  G  CJO,  7’)  such  that  the  solution  v  of  (2.7),  (2.4),  (2.2),  (2.3)  also  satisfies  (2.5). 

PROOF:  As  above,  solve  (2.7),  (2.4),  (2.3)  on  the  triangle  Ao  =  {(f,s)  :  0  < 
M;  a  +  t  <  £°}  and  (2.7),  (2.5),  (2.3)  on  AT  =  {(<,5)  :  0  <  T  -  t,s ;  s  -  t  <  tT  -  T} 
—  obtaining,  in  particular,  v(£°,0)  and  v(T  -  £r, 0).  Then  choose  T  subject  to 
(2.18)  and  specify  Cauchy  data  on  Q  as  above  —  i.e.,  (2.16)  and  v ,  =  0.  Thus,  /  is 
determined  on  Df  —  Pr0  U  P{  U  p{  and  g  on  Pg  =  P90  U  P\  U  P\  with  (2.13)  satisfied 
on  S  =  Pf  n  P9  as  in  the  discussion  above.  Note  that  our  assumptions  give  f',g' 
continuous  on  S. 

Starting  with  the  initial  condition:  £(0)  =  £°,  solve9  (2.12)  to  obtain  I\  This 
solution  curve  satisfies  (2.10)  and  (2.6).  Note  that  the  solution  curve  T  cannot  leave 
S  before  t  —  T  by,  e.g.,  crossing  the  characteristic:  [s  +  <  =  £r  +  7j  —  say,  at  a  point 
(£.,<.)  —  since  this  would  give,  by  (2.10): 

£r  +  T-C  -  £.  =  /(£.  +  /..)  f  g(t.-t.) 

=  /(£r  +  T)  +g{fT  +T  -2t.) 

=  |£T-?(£r-r))  +  9(£r  +  r-2<.), 

contradicting  the  fact  that  |g'|  <  k( 2  <  1/2;  similarly,  T  cannot  leave  S  by  crossing 
|s  -  t  =  £r  -  T\  since  |/'|  <  1/2.  It  follows  that  £(•)  is  well-defined  on  (0, T|  with 
f.(T)  ^  £r. 

This  specification  of  7\  £(■)  defines  Qt ■  In  view  of  (2.6),  the  geometry  of  T  ensures 
that  Qt  C  A0UA7-US  so  /,  q  are  both  defined  on  all  of  Qt-  Now  (2.8)  gives  a  (weak) 
solution  v  of  (2.7)  and  the  construction  ensures  that  we  also  have  (2.4),  (2.5),  (2.2), 
(2.3)  as  desired.  □ 

°A  solution  exists  since  the  light  li.nucl  side  of  (2.12)  is  continuous  in  l  (as  long  as  one  remains  in 
S).  Our  assumptions  do  not  give  a  l.ipschitz  condition;  thus  uniqueness  — ■  not  really  needed  for  the 
Theorem  —  would  have  to  he  obtained  hy  a  different  argument:  noting  that  (2.13)  makes  the  right 
hand  side  of  (2.10)  contractive  in  t. 


Remarks:  The  regularity  required  for  the  data  —  continuity  of  v°,  w°,  vj ,  wT  — 

can  he  weakened  by  a  limit  argument.  In  this  case  one  retains  (2.6)  although  i  may 
no  longer  be  continuous.  There  is  also  no  really  new  difficulty  in  considering  spatially 
variable  coefficients:  pvtt  —  (ftv,)f.  The  argument  must  be  modified  in  detail10  but 
is  essentially  similar  in  spirit. 

ft  is  worth  noting  that  the  construction  in  the  proof  above  can  be  viewed  from  a 
slightly  different  perspective.  After  solving  on  Ao  and  A r  one  fills  in  Cauchy  data  on 
a  ‘RaP  ^  so  as  to  consider  an  initial  value  problem  with  the  roles  of  t  and  s  reversed 
—  taking,  say,  v°,vT  (perhaps  suitably  extended  to  larger  s)  as  boundary  data.  The 
choice  of  data  on  Q  and  the  bounds  (2.19)  are  such  as  to  ensure  that  |t/,|  <  k  <  1 
which  then  makes  the  right  hand  side  of  (2.10)  contractive  in  f.  so  (2.10)  determines 
£(<)  uniquely  for  each  t.  (One  must  also,  of  course,  verify  that  this  gives  6(  )  as  an 
‘admissible  control  function’.) 

Stated  this  way,  it  is  clear  that  the  argument  does  not  really  depend  on  the 
linearity  of  (2.7)  or  on  the  possibility  of  a  decomposition  (2.8).  Without  providing 
details,  we  note  only  that  one  can  treat  (2.1)  from  this  point  of  view  under  quite  mild 
conditions  on  0().  If  we  set  <f>  :=  rp~l,  w  =  vt  and  u  =  0(v,)  so  vf  =  <f>{u)  we  can 
rewrite  (2.1)  in  terms  as  a  system: 


(2.20) 


u,  —  pwt  W,  =  4>{u)t 


suggesting  the  role-reversed  problem.  The  crux  of  the  argument,  now,  is  to  show  that 
‘reasonable  bounds’  on  the  initial  and  target  data  can  be  imposed  to  ensure  that  one 
will  have  (for  a  suitable  construction  of  Q  and  of  the  data  there,  satisfying  m  wt  =  u) 
the  necessary  bound:  |v,|  —  |<£(u)|  <  k  to  be  able  to  apply  the  Contractive  Mapping 
Theorem  to  (2.10).  □ 


3.  TRANSVERSE  VIBRATION 


In  this  section  we  present  an  extremely  brief  introduction  to  the  class  of  problems 
which  might  be  considered.  As  above,  we  consider  the  robot  arm  as  a  nominally 
horizontal  beam  with  a  load  in  at  the  tip  and  with  no  dissipation  or  body  forces  (e.g., 
ignoring  friction,  gravity,  . . . ).  Here,  however,  we  consider  transverse  vibrations.  For 
simplicity  we  consider  only  vertical  vibration,  so  the  principal  unknown  is  the  vertical 
deflection  u  =  n(t,s).  Thus,  the  ‘laboratory  coordinates’  of  a  material  point  s  will 
be  \r.,y,z\  —  [£(f)  -  ,s,0,  u{t ,  *)j.  Furt  her,  we  ignore  shear  and  work  with  a  linearized 


'"Eg.,  the  characteristics  nre  n<>  l.uigci  straight  lines  and  the  propagation  along  them  may  no  longer 
lie  constant.  We  now  must  lake  k  in  !>e  a  constant  less  than  the  minimum  sound  speed  and  impose  a 


more  stringent  condition  (2.10). 


theory.11 

We  begin  with  the  ‘clamped1  boundary  conditions  at  the  prismatic  joint: 

(3.1)  u  =  u,  =  0  at  the  joint:  s  —  £(<). 

Next,  in  this  setting  a  geometric  analysis  shows  that  the  stress  deformation  for  a 
segment  ds  of  the  beam  is  proportional  to  [radius  of  curvature)-2  so  the  strain  energy 
(elastic  potential)  is  ss  *0ti2,ds  where  /?  (sometimes  written  as  El)  corresponds  to 
stiffness.  The  corresponding  kinetic  energy  is: 

fi{t) 

I  /  \p  ( l 2  -I-  u2)  -f  «  u2(j  ds  +  i(m  (£2  -I-  tt2)  +  au2,)|,=n 

•'o 

where  o  is  the  cross-sectional  moment.  By  Hamilton's  Principle,  we  obtain  the  Euler- 
Bernoulli  beam  equation12 

(3.2)  putl  +  (ou,«),  +  (/?«„)„  =  0  on  Qt 
(where  Qt  is  as  above  in  (1.1))  and  the  other  boundary  conditions: 

(3.3)  mutt  -(-  au,„  =  (/?  «,,),,  a  u,tt  4-  /?  u„  =  0  at  the  tip:  s  =  0. 

The  system,  if  the  control  function  f.(-)  were  specified,  would  consist  of  (3.2)  with 
(3.1)  and  (3.3)  together  with  initial  conditions  (2.4).  We  note  here,  however,  that  the 
type  of  problem  having  particular  control-theoretic  significance  is  that  which  arises 
when  £(-)  is  determined  in  feedback  form: 

(3.4)  t(t)  =  A[u(t,)|, 

i.e.,  one  couples  (3.4)  with  the  other  conditions  as  an  implicit  characterization  of  the 
free  boundary  I\ 

In  (3.4),  we  consider  A [ •  j  as  some  suitable  functional  —  presumably  a  function 
u/)  for  some  (finite)  set  of  implementable  observations:13  ut  =  C(u).  For  example, 
one  might  take  A  to  be  of  the  form  A [uj  =  $(C[u|)  with 

(3.5)  C  :  n  {u|,_n,  t o. 

1  (l»  of  a  linearized  stress-strain  relation  is  a  consequence  of  tlie  fundamental  beam  assumption: 
‘infinitesimal’  cross-section,  Riving  infinitesimal  deformations.  The  significant  assumption  is  the  the 
simplifying  but  rather  restrictive  one  that  the  arm  always  remains  straight  enough  to  ignore  the 
so-called  ‘geometric  nonlinearil  ies’  essentially,  that  one  can  adequately  approximate  (radius  of 
curvature)-3  by  u3,,  linear  velocity  by  |fi  t  w,k|,  angular  velocity  |arctanuT|f  by  -uw,  etc. 

l3Note  that  the  equation  is  just  the  standard  beam  equation  when  we  use  a  as  the  independent 
variable.  This  formulation  is  somewhat  different  from  that  of  |2|. 

1 '’Alternatively,  we  might  consider  A  <t>(z)  where  z  is  the  solution  of  an  auxiliary  system  of  the 
form:  z  —  F(z,w). 


Such  an  observation  operator  C  might  be  implementable  by,  say,  visual  observation 
of  the  tip  position  and  placing  an  accelerometer  on  the  beam  at  s*,  allowing  for  the 
(known)  effect  of  the  joint  motion.  This  would  give,  as  an  example, 

(3.6)  f{t)  =  $(u(f,0),u<t(f,s.)) 

for  some  given  function  $. 

Once  a  specific  form  of  the  feedback  functional  A [ •]  has  been  prescribed,  the 
problem  is  analysis14  of  the  full  system:  (3.2),  (2.4),  (3.1),  (3.3),  and,  e.g.,  (3.6)  to 
determine  existence  and  uniqueness  as  well  as  such  properties  of  control-theoretic 
interest  as  stability. 
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